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We analyze several models of atomic detectors in the 
context of the measurement of coherence properties of 
matter waves. In particular, we show that an ionization 
scheme measures normally-ordered correlation functions 
of the Schrodinger field, in analogy with the optical sit- 
uation. However, it exhibits a sensitivity to exchange 
processes that is normally absent in optics. 



Optical coherence theory is based on the observa- 
tion that most quantum measurements that can be 
performed on the electromagnetic field yield a signal 
proportional to normally ordered correlation func- 
tions of that field 0]. A quantized multimode field 
is then said to be coherent to order N if all nor- 
mally ordered correlation functions up to order N 
factorize. No such theory is presently available for 
atomic coherence, probably because until recently it 
had not been necessary to think of atomic samples 
as Schrodinger fields. But the experimental work on 
ultracold atoms, BEC f|-|(| and atom lasers has 
changed that situation, and the need for a proper 
theory of atomic coherence is now quite urgent |pj. 

At least for the case of bosonic fields, it is tempting 
to simply transpose Glauber's coherence theory [[l]. 
This approach has been the de facto situation so far, 
but appealing as it might sound, it must be applied 
with caution, due to the fundamental difference be- 
tween electromagnetic and matter- wave fields. Most 
optical experiments detect light by absorption, i.e. 
by "removing" photons from the light field. This is 
the reason why normally ordered correlation func- 
tions are so important. But atomic detectors work 
in a number of different ways: One can chose to mea- 
sure electronic properties, or center-of-mass proper- 
ties, or both. While one detector might be sensi- 
tive to atomic velocities, another might measure lo- 
cal densities and a third electronic properties only. 
Additional difficulties arise from the fact that atomic 
fields are self-interacting, which significantly compli- 
cates the propagation of atomic coherence as com- 
pared to the case of light. From these remarks, it 
should be clear that a theory of matter waves coher- 
ence is much richer than its optical equivalent. Yet, 
like Glauber's coherence theory, it should be opera- 
tional and based on explicit detection schemes. 

The goal of this note is to analyze several ideal 
atom detectors and to determine which correlation 



functions of the matter-wave field they are sensitive 
to. The systems we explicitly consider are a non- 
resonant atomic imaging system such as used e.g. 
in the MIT BEC experiments, and a detector work- 
ing via atomic ionization. We show that while the 
off-resonance imaging detector is sensitive to density 
correlation functions, a narrow-band ionization de- 
tector measures normally ordered correlation func- 
tions of the Schrodinger field itself, in analogy with 
the optical case. Intermediate situations are more 
complicated, due to the quadratic dispersion of mat- 
ter waves. Higher-order detection schemes also in- 
volve exchange terms usually absent in the optical 
case. 

Nonresonant imaging 

Consider first atomic detection by non-resonant 
imaging, a situation where a strongly detuned elec- 
tromagnetic field interacts with the atoms in the 
sample in such a way that it induces only virtual 
transitions. We consider for concreteness ground 
state atoms described by the Schrodinger field oper- 
ator #(r) with [^(r), #t(r')] = 5(r - r') for bosons, 
and decompose the electromagnetic field into a clas- 
sically populated mode of wave vector ko and polar- 
ization eo and a series of weakly excited sidemodes 
of wave vectors and polarizations q. After adia- 
batic elimination of the upper electronic state of the 
atomic transition under consideration, the interac- 
tion between the Schrodinger field and the radiation 
field is described to lowest order in the side-modes 
by the effective Hamiltonian 

V = % [ d 3 ri^4^* f (r)*(r) 
J °o 

+ ^)% a£e -. t (ko-k,). r \ $t( r )#( r ) (1) 

Oq J 

where k^ is the wave vector of the £-th mode of the 
field, of frequency ue and polarization ee, the sum 
is over all field modes in the quantization volume 
V, and Ei — [huji/ZeoV] 1 / 2 is the "electric field per 
photon" of mode I. The annihilation and creation 
operators at and at satisfy the boson commutation 
relation [o^,oJ,] = S^p. We have also introduced 
the Rabi frequencies fio(r) = d£o(r)(e • eo)/h and 
fit — d£((e ■ ee)/h, and the atom-field detuning 5q = 



1 



uj a — ujq is assumed to be much larger than flo, Sq 
^o(r). 

Assuming that the electromagnetic field is initially 
in the state \£) and the Schrodinger field in the state 
\4> g ) 1 the probability that the system undergoes a 
transition from that to another state is given to first 
order in perturbation theory by 

w= iy^njOi' / d 3 rn (r) [ d 3 r'Q* (r') 
% e/ , J J 

/ dt / dt'&M^Kr'^Mg) 
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where the Schrodinger wave density is defined as 

p(r,t) = *+(r,t)*(r,t) (3) 

and \&(r, t) = U'^(r)U is the time-dependent 
Schrodinger field in the interaction representation 
with respect to the atomic Hamiltonian, i.e. U — 
cxp{—iTLAt/K). 

We further assume for concreteness that all elec- 
tromagnetic sidemodes are initially in a vacuum. 
The measurement on the Schrodinger field is then 
carried out by detecting photons scattered by the 
atoms into the sidemodes, in a fashion familiar from 
resonance fluorescence experiments. The most im- 
portant non-trivial contribution to the fluorescence 
signal is proportional to the intensity |i7o| 2 of the 
incident field, 

w= ^j2\n i \ 2 J d 3 r J d 3 r' dt dt' 

e i(k -k e )-(r-r') e -i(uj -oJi)(t-t') 

x (^|p(r,i)p(r',t')|^), (4) 

and hence is sensitive to the second-order correlation 
function of the sample density. This is to be com- 
pared to the results of Javanainen M , who showed 
that the the spectrum of the scattered radiation is a 
function of (p(r, 0)/5(r, t)}. Indeed, it can be shown 
in all generality that any measurement involving the 
electromagnetic field scattered by the atomic sample 
under conditions of off-resonant imaging are deter- 
mined by correlation functions of the Schrodinger 
field density. 

Ionization 

The reason off-resonant imaging yields a signal de- 
pendent on p(r, t) is that the electric dipole interac- 
tion is bilinear in the Schrodinger field operators. 
This difficulty can however be eliminated if, instead 
of making measurements on the radiation field, one 
detects the atoms directly. One scheme that achieves 



this goal is the ionization method that we now dis- 
cuss. 

Consider a detector consisting of a tightly focussed 
laser beam that can ionize atoms by inducing transi- 
tions from their ground electronic level \g) to a con- 
tinuum level [J The corresponding single-particle 
Hamiltonian is 

H = H cm + H el + V(r) =H + V (5) 

where H cm is the center-of-mass Hamiltonian, H e i 
the electronic Hamiltonian, and V(r) describes the 
electric dipole interaction between the atom and the 
ionizing laser field. H e i has eigenstates <p n and eigen- 
frequencies u> n , H e i\ip n ) = hu> n \ip n ). The corre- 
sponding atomic manybody Hamiltonian is 

Ho = J d 3 r&(r)H V(r) (6) 

where in the Born-Oppenheimer approximation 
^(r) is a multicomponent field with components 
*«(r). 

We are interested in measuring properties of the 
ground state component Vf ff (r) of this field, which 
is dipole-coupled to continuum states ^i(r). We 
assume for simplicity that the center-of-mass wave 
function of these latter states is well described by 
plane waves of momentum q, so that TL may be ex- 
pressed as 

Ho = H g + y^%, (7) 

i 

where 

H i = hJ2(u q + Ui)b\ ti b (l ,i- ( 8 ) 
q 

Here we expanded ^(r) in plane waves as ^E'i(r) = 

Eq^tqW^qi witn b q' = ^qq" 5 "' ) aild = 

Tiq 2 /2M. (Note that the inclusion of ground state 
collisions is straightforward and does not affect our 
conclusions.) 

In terms of the components ^ n ( r ) of the 
Schrodinger field, the electric dipole interaction 
Hamiltonian is 

V = hJ2 J d 3 rfli(r)^l{r)^ g (r)+H.c., (9) 

i 

where fii is the Rabi frequency between the levels 
\g) and and the laser field is treated classically. 



Hot wire detectors can be modeled in a similar fashion. 
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In this detection scheme, one extracts information 
about the state of the field & g (r,t) by measuring, 
e.g. the number of atoms in the continuum. For 
atoms cooled well below the recoil temperature and 
tightly focused laser beams, the spatial size of the 
atomic wave function is much larger than the laser 
spot and we can approximate the electric field E(r) 
by E(r) ~ E5(r - r ), so that Eq. (||) becomes 

V = hJ2^i( r o)M(r )^ g (r ) + B.c. (10) 

i 

We take the atomic system to be initially in the 
state 

|^(o)) = |{VWo)},^(o)). (ii) 

To first order in perturbation theory, the transition 
probability away from that state during the time 
interval At is 

pAt pAt 

w = V|Oi(r )| 2 / dt dt' 
i, q Jo Jo 

<VWO)l**(ro,#!(r ,OlVWO)> 
x (^(0)|4t(r ,i)* 9 (r ,t')|V' 9 (0)} + c.c. (12) 

There is a fundamental distinction between the 
present situation and Glauber's photodetection the- 
ory, because in the present case both the detected 
and detector fields consist of matter waves. There 
is a complete symmetry between these two fields so 
far, and their roles are interchangeable. In order to 
break this symmetry and to truly construct a de- 
tector, we now make a series of assumptions on the 
state of the detector fields ^(r, t). Physically, this 
amounts to making a statement about the way the 
detector is prepared prior to a measurement. Specif- 
ically, we assume that all atoms are in the ground 
state, , i'i(r )|^ ijq (0)) = |0), and that any atom in 
an ionized state will be removed from the sample in- 
stantaneously. In that case, the second term in Eq. 
( |l2| ) vanishes and we have 

pAt pAt 

w = y2\n i (r )\ 2 dt dt' 
i Jo Jo 

^^-''^(roKM 

x^-t'H^^Mifgira,*)^). (13) 

At this point, it is convenient to distinguish three 
different operating regimes: In the first one, only one 
final electronic state is considered, and in addition 
a velocity selector is used to filter just those ionized 
atoms with a given center-of-mass momentum. We 
call this a narrowband single-state detector. The sec- 
ond scheme allows for a broader velocity filter, but 



still considers a single continuum electronic state, 
and we call it a broadband single-state detector. Fi- 
nally, we also discuss a general broadband detector 
where neither the final momentum state nor the fi- 
nal electronic state is narrowly selected. 

More precisely, a narrowband single-state detector 
includes a velocity selector with a bandwidth Aq 
around a central value qo such that for the detection 
times At of interest, one has AtAutq <C 1, where 
Aw q = ?iqoAq/2Af. In this case and for a stationary 
Schrodinger fields Eq. ( |l3| ) reduces to 

Aw 3 

r n& (w,w qo ) = — ^\n(r )\ 2 

pAt 

x / dre-^+^G^-r; r ,r ), (14) 
Jo 

where we dropped the index i of the observed contin- 
uum state for clarity, introduced the ionization rate 
r„(,(w, w q ) = w n b(uj,ojq)/ At and defined the atomic 
normally ordered first-order ground state correlation 
function 

G A (t,t';r ,v ) = (0 fl |*+(ro,i)* ff (ro,Ol0fl)- 

From the Wiener-Khintchinc theorem, we recognize 
that for large enough At, the detector measures the 
spectrum of the Schrodinger field ^ g (ro, 0). 

In the case of a broad single-state detector, in con- 
trast, we have 

pAt 

r lb ~ \n(r Q )\ 2 / dTe-^ T G pr (0,T;r Q ,r ) 
Jo 

x G A (0,T;r Q ,r ) (15) 

where we have introduced the center-of-mass prop- 
agator 

Gp r (*i,i2;ri,r 2 ) = Y / M^W^)e^- tl) . 
q 

(16) 

In that case, the ionization rate is propor- 
tional to the Fourier transform of the product 
of G pr (0, r; ro, ro) and the correlation function 
G^O, r; ro, i"o), or in other words to the convolution 
of the Fourier transforms of these functions. The 
Fourier transform of the center-of-mass propagator 
can therefore be interpreted as the spectral resolu- 
tion of the detector. 

We finally turn to the case of a general broad- 
band detector, where the spectrum of the detector 
is much broader than the spectrum of the detected 
quantity. Assuming that the spectrum of the atomic 
correlation function is centered at a), we find 

r bb ~J7(ro)G A (0,0;ro,ro), (17) 
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where we have introduced the "detector efficiency" 

*?(r )= / dr^|a i (r )| 2 (v|/ l (ro,r)^(r ,0)}e-^. 

i 

(18) 

As expected, a broadband detector is not able to re- 
solve any spectral feature of the Schrodinger field, 
and only measures the local atomic density. 

Higher-order correlations 

The detection of higher-order correlation func- 
tions of the Schrodinger field can be achieved by a 
straightforward generalization of the ionization de- 
tector. For instance, second-order coherence mea- 
surements can be carried out by focussing the laser 
at two locations ri and r 2 , in which case 

V = ^ £fi i (r / ,)$J(r /t )$ fl (r,,)+/r.c. (19) 

(11=1,2 t 

Assuming as before that the continuum states are 
initially empty and for a general broadband detector, 
the joint probability to ionize an atom at ri and 
another one at r 2 is then 

pAt pAt 

w 2 ^ 77(1-1, r 2 )77(r 2 ,ri) / dt x / dt 2 
Jo Jo 

(*J( r i'*i)*l( r 2'*2)* g (r 2! ti)* g (ri,t 2 )) 

pAt pAt 

+ ry(ri)?7(r 2 ) / dh / dt 2 
Jo Jo 

(¥ g (r 1 ,t 1 )i>l(r 2 ,t 2 )i' g (r 2 ,t 2 )i> g (v 1 ,t 1 )), (20) 

where we have introduced the detector cross effi- 
ciency rj(ri,r 2 ) as a straightforward generalization 
of Eq. (fig). The first term in Eq. ( J20j ) is an ex- 
change term resulting from the fact that the detector 
field is a single Schrodinger field. It results from the 
interference between detectors at points ri and r 2 . 
The second term is the usual term also appearing 
in the double photo-detection of optical fields. In 
that latter case, the exchange term does not appear 
because the two detectors used to measure the field 
are taken to be distinguishable. Note also that in the 
position measurement scheme proposed in Ref. [ flO|| , 
interferences do not occur as the set of states ionized 
at each location are taken to be distinguishable. We 
finally remark that as a consequence of the exchange 
term, the signal cannot simply be expressed in terms 
of correlation functions of p(r, t). 

In summary, we have analyzed several detectors 
that permit to access different classes of correla- 
tion functions of the Schrodinger field. Most in- 
teresting perhaps is the ionization scheme, which is 
closely related to the detectors familiar from the de- 
tection of optical fields. However, it presents new 



features, and is in particular sensitive to exchange 
processes. But ionization detectors make destructive 
measurements. This is in contrast to off-resonant 
imaging, which is nondestructive but measures den- 
sity correlation functions instead of the more fa- 
miliar normally-ordered correlation functions of the 
Schrodinger field itself. 
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